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Abstract 

Random walk models in one-dimensional disordered media with an oscillatory input current are 
investigated theoretically as generic models of the boundary perturbation experiment. It is shown 
that the complex admittance obtained in the experiment is not self- averaging when the jump rates 
Wi are random variables with the power-law distribution ~ Wi'^^^ (0 < a < 1). More 

precisely, the frequency-dependence of the disorder-averaged admittance (x) disagrees with that of 
the admittance x of any sample. It implies that the Cole-Cole plot of (x) shows a different shape 
from that of the Cole-Cole plots of x of each sample. The condition for absence of self-averaging is 
investigated with a toy model in terms of the extended central limit theorem. Higher dimensional 
media are also investigated and it is shown that the complex admittance for two-dimensional or 
three-dimensional media is also non-self-averaging. 

PACS numbers: 02.50.Ey, 05.40.Fb, 71.23.Cq, 72.80.Ng 
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I. INTRODUCTION 



Disordered systems like amorphous semiconductors have a large number of remarkable 
properties inherent in their disorder, like the dispersive transport Hence, the disor- 
dered systems have wide areas of application like the photoreceptors of solar batteries and 
photocopiers^. So, they attract continuous research interests. In investigation of such re- 
markable properties, two kinds of statistics are to be considered. One is of course statistics 
of thermal fluctuation, which is considered in traditional statistical physics. The other is 
statistics of random structure of disordered systems. In order to treat this type of statis- 
tics, random structure is described with random variables and the probability distribution 
function of the random variables characterizes the substances. 

It is usually assumed in studies of disordered systems that a sample used in experiments 
is sufficiently large. Hence, a measurement of any observable in such a system corresponds 
to an average over the ensemble of all realizations of the disorder. Quantities for which this 
assumption is valid are said to be self-averaging. The assumption means no sample depen- 
dence in the quantities. So, if the assumption is valid, reproducibility of experimental results 
for different samples is guaranteed. In addition, if the assumption is valid, in theoretical 
analyses only a disorder-averaged quantity is to be calculated since it should coincide with 
experimental results for any sample. 

The expectation that the assumption of self-averaging is valid is based on the law of large 
numbers. The theorem says that the mean of independent random variables is equal to the 
expectation value with probability one. However, disordered systems such as amorphous 
semiconductors are strongly disordered, i.e., the expectation value of the random variables 
that describe structural disorder diverge. For example, the dispersive transport has been 
successfully explained by strongly disordered random variables- . In the theory of the disper- 
sive transport, the transport is modeled by a random walk model on the disordered lattice 
where the dwell times of carriers at lattice sites are strongly disordered random variables. 

The distribution of strongly disordered random variables is broad. Hence, properties of 
limited areas of samples like the maximal dwell time govern the macroscopic properties like 
the macroscopic mobility. One simple example is the finite contribution of the maximal 
term of a set of strong disordered random variables to the sum of the set (see Eq. fl81|) 
in appendix.). The properties of the limited areas with principal contribution fiuctuate 
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from one sample to another. Hence, the macroscopic properties of such materials may 
fluctuate. Thus, the assumption of self-averaging should be tested for its validity in strongly 
disordered systems. In fact, there are several reports on absence of self-averaging in transport 
phenomena in disordered media: Ref. 3 for mean square displacement, ref. 4 for the mean 
first passage time in the Sinai model. 

In this paper, we study theoretically self- averaging properties of physical quantities ob- 
tained in a boundary perturbation experiment. The experimental technique was introduced 
recently to measure the optoelectrical properties of amorphous semiconductorsi^. In this 
experimental method, an oscillatory perturbation with laser light is applied at one end of 
a sample and a response of the photocurrent from the other end is measured. We use as 
a generic model of the boundary perturbation experiment a random walk model in a one 
dimensional lattice with a boundary condition oscillating in time. The potential surface 
that carriers in the amorphous semiconductors feel is rugged due to structural disorder and 
it is called a rugged energy landscape. Consequently, the jump rates of the carriers are ran- 
dom variables and its probability distribution characterizes the disorder. It is known that 
the random jump rates obeys the power-law distribution with a negative exponent^ and it 
implies that the energy landscapes of the disordered media are extremely rugged. 

It is shown that the complex admittance obtained in the experiment is not self-averaging 
when the jump rates Wi are random variables distributed by the power-law distribution with 
the negative exponent. More precisely, the frequency-dependence of the disorder- averaged 
admittance (x) disagrees with that of the admittance x of ^iny sample. It implies that the 
Cole-Cole plot of (x) shows a different shape from that of the Cole-Cole plots of x of sach 
sample. The condition for absence of self- averaging is investigated with a toy version of the 
complex admittance. The complex admittance for higher dimensional media is investigated 
and shown to be non-self-averaging. 

The rest of this paper is organized in the following way: In Sec. m the random walk model 
used in the present study is explained. In terms of random walk, the complex admittance 
is related to a more physically relevant quantity to describe the transport in a disordered 
medium in Sec. IIIII In Sec. IIVI absence of self-averaging for the complex admittance of one- 
dimensional media is shown. In Sec. |3 absence of self-averaging is displayed in the Cole- 
Cole plot impressibly. In order to understand the condition for absence of self-averaging, 
a toy model is studied with the extended central limit theorem in Sec. IVII Absence of 
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self- averaging for higher dimensional disordered media is shown in the similar way for the 
one-dimensional case in Sec. I VI II Discussions and conclusions including consideration of 
possibilities that non-self- averaging properties are observed in real experiments are given in 
Sec. IVlllI A part of this work has been reported in the brief note^ without consideration 
based on the probability distribution of the normalized mean first passage time, study of 
the toy model and the results for higher dimensional media. 

II. A GENERIC MODEL OF THE BOUNDARY PERTURBATION EXPERI- 
MENT 

Response of a system to an external field is a standard tool to be utilized in the study of 
condensed matter. Recently, the intensity modulated photocurrent spectroscopy^ has been 
introduced to investigate transport properties in amorphous semiconductors. In the experi- 
ment, the electrode is illuminated with a laser and the incident light intensity is harmonically 
modulated at frequency u. The light absorbed by an optical transition generates the pho- 
tocurrent. The light intensity absorbed in the electrode consists of a background intensity 
and an oscillating component with small amplitude (pioj)., which respectively gives rise to a 
steady state photocurrent and a harmonically varying photocurrent /\i{uj). Since the har- 
monically varying photocurrent Ai(u;) may show a phase shift with respect to the absorbed 
light flux the optoelectrical admittance /\i{uj) / becomes a complex number. 

Since in this experiment an oscillatory perturbation is applied at one end of a system and 
a response is measured, the experimental technique belongs to a generic method called the 
boundary perturbation method (BPM). In the presence of a periodically forced boundary 
condition on one end of the system, the output from the other end of the system is in pro- 
portion to the perturbation in the linear response regime and the proportionality constant is 
called the admittance. The frequency dependence of the admittance contains various infor- 
mation of the dynamics of the system. The BPM is expected to provide useful information 
on transport properties of disordered media. 

We introduce as a generic model of the BPM a random walk model in one-dimensional 
disordered lattice segment of + 1 sites with an oscillatory input current (see Fig. Q). The 
lattice sites are denoted by integers, n = 0, 1, . . . , A^. The dynamics of a random walking 
particle can be described by a master equation for the probability Pn{t) that the particle is 
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FIG. 1: The structure of the one-dimensional lattice segment. The oscillatory current perturbation 
J{t) is introduced at the site 0. The site is assumed to be a sink. 

at the site n at time t > 0. The master equation for the model is written as 
dPnit) 

if) - {Wn-1 n + Wn+l n)Pn{t) + Wn n+lPn+lif) , (l) 

where Wm n denotes the random jump rate of a particle from the site n to the site m. Here, 
we assume that the random walking particle jumps only to the nearest neighbor sites. The 
probability distribution of jump rates characterizes the disordered medium. We introduce a 
perturbation at the left end, so that the equation for Po{t) is given by 

= -w^ oPo(t) + wo iPi(t) + J{t). (2) 

J{t) is the oscillatory current perturbation at the site 0. The right end of the system is 
assumed to be a sink and the equation for Pn-i is given by 

^ = Wn-1 N^2PN-2{t) - {wn-2 N-1 + Wn N-l)PN~l{t). (3) 

Since the current perturbation J{t) oscillates in time around a positive average with the 
amplitude A J, the response of the output current wn N-iPN-i(t) oscillates around its 
stationary state with the amplitude wn n-i^Pn-i at the same frequency with a phase- 
shift. The admittance is defined by the ratio of these two amplitudes 

/ N Wn n-i'^Pn-1 
Xn{uj) = — . (4) 

The system of equations to determine the amplitudes AP„ of oscillation of the probabilities 
is derived from the master equation Eqs. (HJ), Q and 0; 

iujAPn = Wn n-l'^Pn-l " {Wn-1 n + Wn+1 n)APn 



+Wn n+1 AP for n = 2, 3, . . . , AT - 2, 
ztuAPo = -wi oAPo + lAPi + AJ, 
iujAPn-i = wn-1 n~2^Pn~2 - {wn-2 n-1 + wn n~i)^Pn^i- (5) 

III. THE COMPLEX ADMITTANCE AND THE FIRST PASSAGE TIME DIS- 
TRIBUTION 

We show below that the admittance can be related to the first passage time distribution 
Pn o{t) (FPTD), which is the probability density that the particle which starts at the site 
at time arrives for the first time at the site at time t. Since the site is a sink, the 
FPTD is given by the output current from the site when there is no input current and 
the particle starts from the site at time 0, i.e. 

Fn oit) = wn N^iPw-iit). (6) 

The Fourier-Laplace transform of the master equation for the case is written as 

iuJpn{iuj) = Wn „_lP„_l(ia;) - [Wn-l n + Wn+1 n)Pn{i^^) 

+Wn n+iP{i^) for n = 2, 3, . . . , A^ - 2, 

iujPQ{iuj) = -wi oPiiuj)o + Wo iPi{iuj) + 1, 

iujPN^l{iuj) = Wn~1 N-2PN-2{ii^) - {wn~2 N~l + 

Wn N-i)PN-i{iuj). (7) 

This set of equations is identical to Eq. (0) divided by A J. It means that AP„/AJ = 
Pn{iuj) and hence the Fourier-Laplace transform of the FPTD q{s = iu) is identical to 
Wn Ar-iAP/v„i/A J. Thus, from Eqs. (jU and (0), we conclude that the admittance is equal 
to the Fourier-Laplace transform of the FPTD; 

XAr(u;) = Fn oii^^)- (8) 

We make the low-frequency expansion of the admittance to see the behavior near the 
static limit. Since the admittance is given by the Fourier-Laplace transform of the FPTD, 
the admittance at zero frequency equals to unity due to normalization of the FPTD; 

Fn o(0) = / dtFN o{t) = 1- (9) 
Jo 
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The mean first-passage time (MFPT) tjv, which is the first moment of the FPTD, is given 
by = « dx{uj = 0)/duj. Thus, the low-frequency expansion of the admittance is given by 

XNiuj) = l~itN UJ + 0{uj'^). (10) 

The MFPT is given by^ 

N-l ^ N-2 ^ N~l i 

^-=E^+E^E n (11) 

k=0 ^k-l k k=0 k i=k+l j=k+l j 

IV. ABSENCE OF SELF- AVERAGING FOR THE COMPLEX ADMITTANCE 

Since we are interested in the complex admittance for amorphous semiconductors, we 
employ the power law probability distribution of the random jump rates Eq. (fT^ . which is 
valid for amorphous semiconductors^: 

aw°'~^ if < w < 1 

(12) 

otherwise. 

There is possibility of the absence of self-averaging when disorder of the random dwell 
time at a site is strong as suggested in Sec. HI Since the dwell time is proportional to inverse 
of the jump rate, we are interested in the case where the first inverse moment of the random 
jump rate (l/w) diverges, i.e. < a < 1. It is important to note that the inequality 
< a < 1 holds for amorphous semiconductors^. 




A. The probability distribution of the MFPT 

In this subsection, we analyze the probability distribution of the MFPT Eq. (jllj) . Since 
the MFPT gives the coefficient of the low frequency expansion of the admittance, the results 
of the analysis on the MFPT make it possible to determine whether the admittance is 
distributed with finite variance, i.e. non-self-averaging, or not. We show below that the 
complex admittance Eq. (llUj) is non-self- averaging when < a < 1. 

In order to analyze Eq. (jlip . we employ the random trap model^^, where the jump rates 
Wij depends only on j. For the random barrier model, where the jump rates have the 
symmetry such that Wij = Wji, the quantitatively same results for self-averaging properties 
are obtained. We discuss below the probability distribution function of the normalized 
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The normalized MFPT 

FIG. 2: The probability distribution function of the normahzed MFPT /N^^'^'^^ for the random 
trap model where the probability distribution function of the random jump rates is power-law Eq. 

with exponent a = 1/2. The distribution is computed from 500000 samples of disordered 
chains of the length = 2, 5, 10, 100, 300, 1000. It is clearly seen that the probability distribution 
functions lie on a same curve when the length of lattice is sufficiently large. 

MFPT for the three cases of the value of a. It is important to note that the three cases 
exhaust all possibilities of the value of a. The choice of the normalization constants is 
suggested by the extended central limit theoren>i2,. The theorem is summarized in appendix 
□ 

1. The case when < a < 1 is considered. It is numerically shown that the probability 
distribution function of the normalized MFPT defined as tjq /N^/°'^^ converges to a dis- 
tribution function with non-zero variance in the limit ^ oo. Namely the normalized 
MFPT, limTv^oo ^Tv/^^^"^^, is a random variable even for the infinitely long chain. In 
addition, the distribution function behaves asymptotically as p(x) ~ x~'^{l < 7 < 2) 
for large x. It is shown in Sec. IVllBl that 7 = a + 1. Figure |21 shows the probability 
distribution function of ij\f/N^^"^^ when a = 1/2 in Eq. (fT^ . It is clearly seen that 
the probability distribution functions for sufficiently long chains lie on a same curve. 
Figure El shows the tail of the probability distribution which appears in Fig. |21 It 
behaves asymptotically as p{x) ~ x~'^ (7 = 1.5 < 2). 
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The normalized MFPT 



FIG. 3: The asymptotic behavior of the probabiUty distribution function of the normahzed MFPT 
computed from 5000000 samples. The soHd hne denotes the distribution function when N = 1000 
and a = 1/2 and the dashed hne represents the function proportional to x~^'^. One sees clearly 
that the distribution function behaves asymptotically as p{x) ~ x~'^{'y = 1.5 < 2). 

2. When a = 1, it is numerically shown that the probability distribution function of the 
normalized MFPT defined as In/N^ - &7v/2 converges to a distribution function with 
non-zero variance in the limit ^ oo. Here, the centering constant is defined as 



In appendix El we show that diverges as — *■ oo. Figure. |3] shows the probability 
distribution function of the normalized MFPT when a = 1. It is clearly seen that the 
probability distribution functions lie on a same curve when the length of lattice is 
sufficiently large. 

3. When 1 < a, it is also shown numerically that the probability distribution function 
of the normalized MFPT defined as t^/N'^ converges to the Dirac delta function with 
the support at {l/w)/2. FigureElshows the probability distribution function of i]\f/N^ 
when a = 1.5. One sees that the probability distribution concentrates on the point 
{l/w)/2 = 1.5 when the length is large. 




(13) 
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The normalized MFPT 

FIG. 4: The probability distribution function of the normahzed MFPT t]\i/N'^ — /2 for the ran- 
dom trap model where the probability distribution function of the random jump rates is power-law 
Eq. (|12p with exponent a = 1. The distribution is computed from 500000 samples of disordered 
chains of the length = 2, 5, 10, 100, 300, 1000, 3000. It is clearly seen that the probability distri- 
bution functions lie on a same curve when the length of lattice is sufficiently large. 

B. Absence of self-averaging in the low frequency behavior of the complex admit- 
tance 

By using the results for the probability distribution of the normalized MFPT, we discuss 
below frequency dependence of the complex admittance Eq. (fTU|) for the three cases of the 
value of a: < a < 1, a = 1 and 1 < a. 

1. When < a < 1, it is shown that the complex admittance is non-self-averaging. In 
order to see frequency dependence of the admittance, the long chain limit (A^ oo) 
should be taken with keeping i^oj finite. Otherwise, the admittance becomes a trivial 
constant or 1. Since it was shown that the normalized MFPT ijq/N^^"^^ is finite 
with probability one, the long chain limit should be taken with fixing ojN^/"'^^ finite. 
Thus, from Eq. the admittance for the infinitely long chain is given by 

lim Xn{i^) = 1 — i t + . . . (14) 
10 




FIG. 5: The probability distribution function of the normahzed MFPT t^/N'^ when a = 1.5. It is 
computed from 500000 samples of disordered chains of length N = 100, 300, 1000, 10000. It is seen 
clearly that the probability distribution concentrates on the point {l/w)/2 = 1.5 when the length 
of lattice is sufficiently large. 

where i is the normalized MFPT defined as t = liniTv^oo ^Af/^^^"^^ and lim jv^oo 
means the long chain limit with fixing luN^^""^^ at the finite value u. Equation (jl4|) in- 
dicates that the admittance for the infinitely long chain is a function of the normalized 
MFPT if: and the scaled frequency u. From our result for the probability distribution of 
the normalized MFPT when 0<a<l,tisa random variable. Thus, the admittance 
is also a random variable and hence the admittance is non-self-averaging. 

2. When a = 1, it is shown below that the admittance is self- averaging in the low 
frequency range. By introducing the normalized MFPT t^- defined as 

i^uj is rewritten as 

2 

tjvuj = N%NUj{—iN + l)/2. (16) 
On 

It was shown in Fig. |3]that the normalized MFPT tN is finite with probability one. 
Thus, in order to see frequency dependence, the long chain limit should be taken with 
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keeping the scaled frequency cD = finite. Since diverges as ^ oo, the 

complex admittance for the infinitely long chain is given by 

2 ~ 

lim xn{uj) = lim 1 - iN^nuj l2i—tN + I) + ■ ■ ■ 

N^oo iV— >oo Oat 

= l-iCd + ... (17) 

Since Eq. (jl7|) has no sample- dependence, the complex admittance is self- averaging in 
the low frequency region. 

3. When 1 < a, it is shown below that the complex admittance is self-averaging. In 
this case, it was shown that the scaled MFPT tN/N"^ is equal to {l/w)/2 in the limit 
N oo. Thus, the admittance for the infinitely long chain is given by 

lim XNico) = l-l-(-)u + ... (18) 



N^oc 2 \W 

ui—uiN-' 



It means that the complex admittance is self-averaging. 
From Eq. (|14j). the low- frequency expansion of the disorder-averaged admittance is given 



by 



]im xn{oo)) = l-t{t)Co + ... (19) 



When < a < 1, the distribution function of the normalized MFPT t behaves asymptot- 
ically as p{x) ~ x~'^. It is shown later in Sec. IVIIBI that 7 = a -|- 1 and it means that 
1 < 7 < 2. Hence, the expectation value (t) of the normalized MFPT diverges. Since it im- 
plies that the coefficient of the first order of low u expansion diverges, the disorder-averaged 
admittance is a non-analytic function of iuj which must behave as 

( lim xn{uj)) - 1 ~ {tCoT (20) 

where < ;U < 1. 

In order to test the foregoing observation, we numerically solved Eq. (0) for one- 
dimensional chain of 20 sites. The jump rate Wi j depends only on j (the random trap model) 
and the jump rates are generated from the power-law distribution Eq. (jl2j) . Figure El shows 
the low-frequency behavior of the real part and the imaginary part of the admittances when 
a = 1/2. One clearly sees that the low-frequency behavior of the disorder-averaged admit- 
tance is completely different from the admittance for any of the samples. On the other hand, 
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The scaled frequency The scaled frequency 

FIG. 6: The frequency dependence of the real part (a) and the imaginary part (b) of the complex 
admittances when the exponent of the probability distribution of jump rates is a = 1/2 and the 
length of a chain is = 20. Solid lines represent the admittances for each of 100 samples. The 
dashed line represents the admittance averaged over 5000 samples. Though the admittance of each 
sample is proportional to the scaled frequency a), the averaged admittance is proportional to 
(/" ^ 1/2). 
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FIG. 7: The frequency dependence of the real part (a) and the imaginary part (b) of the complex 
admittances when a = 1.0 and the length of a chain is 20. Solid lines represent the admittances 
for each of 100 samples. The dashed line represents the admittance averaged over 5000 samples. 
Both the admittance of each sample and the averaged admittance are proportional to the scaled 
frequency Co. 
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FIG. 8: The frequency dependence of the real part (a) and the imaginary part (b) of the complex 
admittances when a = 1.5 and the length of a chain is 20. Solid lines represent the admittances 
for each of 100 samples. The dashed line represents the admittance averaged over 5000 samples. 
Both the admittance of each sample and the averaged admittance are proportional to the scaled 
frequency a). 

figures IZI and IHl show the low-frequency behavior for a = 1.0 and for a = 1.5. One clearly 
sees that the admittances for individual samples and their disorder-averaged admittance are 
both proportional to u. Although sample- dependence is observed, the sample-dependence 
is considered as the finite-size effect or effects from the terms of 0{uj'^). 

V. NON-SELF-AVERAGING IN THE COLE-COLE PLOT 

In the literature of experiments of the BPM^, the Cole-Cole plot of the admittance is 
employed for the analysis of experimental results. It is a parametric plot of the imaginary 
part of the admittance against the real part. Figure IHl shows the Cole-Cole plot of x and 
(x) obtained by solving Eq. ^ when a = 1/2 in Eq. (fT^ . One sees that the shape of the 
Cole-Cole plot of the disorder-averaged admittance is completely different from the shape 
of the Cole-Cole plot for each of samples. The boundary of the region where the Cole-Cole 
plot of the admittance of each sample scatters is the Debye semi-circle, whose center is at 
(1/2,0) and the radius is 1/2. In Fig. El solid circles denote the values of the admittance 
at a given frequency and a diamond is the average over these values. Since the admittance 
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FIG. 9: The Cole-Cole plots of the admittances when the exponent of the probability distribution of 
jump rates is a = 1/2 and N = 100. Solid lines represent the admittances for each of 30 samples. 
The dashed line represents the admittance averaged over 5000 samples. The red solid circles 
represent the admittance of each sample at tD = 0.5. Although the red solid circles scatter outside 
the Debye semi-circle, the yellow diamond which denotes the averaged admittance at a) = 0.5 is 
inside of the Debye semi-circle. 

of each of samples at the same frequency (a solid circle) scatters around the Debye semi- 
circle, the admittance averaged at the same frequency (the diamond) comes inside the Debye 
semi-circle. 

In order to analyze this non-self-averaging property of the admittance in the Cole-Cole 
plot, we prove rigorously that the Cole-Cole plots of each samples must lie outside the Debye 
semi-circle. We first note that the inverse of the admittance satisfies the following recursive 
relation: 

Xn (^) = + 1 + Xn-l(^) X„-2(^), (21) 

where x~^(u;) is the inverse of the admittance for a chain of length n. The recursive Eq. 
(j21|) is derived from the following two equations for the FPTD obtained from its definition: 
-^.+1 i{s) = ipi+i i{s)+ijji-i i{s)Fi+i i_i(s) and F^+i i_i(s) = F^+i i{s)Fi i_i(s), where Fi j{s) 
denotes the Laplace transform of the FPTD from site j to site i and ipi denotes the 
Laplace transform of the waiting time distribution for the jump from site j to site i. The 
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recursive relation proves inductively that the inverse of the admittance obeys the following 
four inequalities in the frequency region [OjCiJn], where Un is the positive smallest zero of the 
real part of and is equal to the positive smallest zero of Xn'- 

x-J" > 

"»-■' " " (22) 

Xn ^ Xn— 1 

-1' -1 ' 

Xn ^ Xn—1 • 

Since the second inequality shows c<j„ < ujn~i, the fourth inequality imphes Xn^' < Xn-i < 
■■ ■ < Xi = 1 in the region [0,t^„]. Since Xn = Xn'/{Xn + Xn" ), the inequality Xn < 1 
means {xn — 1/2)^ + Xn"^ > (1/2)^- Thus, the Cole-Cole plot of the admittance of each 
sample can exist only outside the Debye semi-circle when the frequency is smaller than the 
positive smallest zero of Xn- 

On the other hand, we can show in the following way that the Cole-Cole plot of the 
averaged admittance is located inside the Debye semi-circle when < a < 1. We consider 
the angle 6 between the tangent of the Cole-Cole plot at (1, 0) and the horizontal axis. From 
Eq. (j^Uj) . 6 = /U7r/2, which is less than tt/2 since /i < 1, is obtained. However the angle 6 
for the Debye semi-circle is equal to 7r/2. It means that the Cole-Cole plot of the averaged 
admittance exists inside the Debye semi-circle when < a < 1. Since the Cole-Cole plot of 
the admittance of any sample exists outside the Debye semi-circle, the Cole-Cole plot clearly 
shows non-self-averaging property of the admittance. 

When a > 1, the angle 6 of the Cole-Cole plot of the averaged admittance is 7r/2 and 
is equal to that of the Cole-Cole plot of each of samples as seen in Figs. ^1 and ^2 This 
property of the angle 6 is derived from the regular behavior as (x) — 1 ~ icD (see Eq. (II 7|) 
and Eq. ^ .). 

It is important to note that Fig. 021 suggests that the complex admittance is non-self- 
averaging when a = 1. One sees that the disorder- averaged complex admittance (x) exists 
at different position from that of the complex admittance for each samples. Although it 
was shown that the low frequency behavior is self-averaging, it seems that the complex 
admittance at the finite frequency is non-self-averaging. Consequently, we conclude that the 
complex admittance is non-self-averaging when < a < 1. 
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FIG. 10: The Cole-Cole plots of the admittances when the exponent of the probability distribution 
of jump rates is a = 1.0 and N = 100. Solid lines represent the admittances for each of 30 
samples. The dashed line represents the admittance averaged over 5000 samples. The red solid 
circles represent the admittance of each sample at a) = 0.68. The green diamond which denotes 
the averaged admittance at lD = 0.68 exists at different position from that of the red solid circles. 

VI. UNDERSTANDING OF NON-SELF- AVERAGING OF THE COMPLEX AD- 
MITTANCE WITH A TOY MODEL 

In the previous section, we discussed absence of self-averaging for the admittance when 
disorder of the jump rates is strong. In this section, in order to consider the conditions for 
absence of self-averaging we introduce a simple model, which is a toy version of the complex 
admittance and can be analyzed rigorously. This simple model will provide understanding 
of the reason why the admittance is non-self-averaging, since absence of self-averaging for 
the simple model is shown below with the same discussion as employed in the analysis of 
the complex admittance. 

Analysis of the admittance is difficult because it does not consist of a mean of indepen- 
dent and identically distributed random variables. Hence, we consider a toy model of the 
admittance which is a simple function of a mean of independent and identically distributed 
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FIG. 11: The Cole-Cole plots of the admittances when a = 1.5 and N = 100. Solid lines represent 
the admittances for each of 30 samples. The yellow dashed line represents the admittance averaged 
over 5000 samples. The red solid circles represent the admittance of each sample at a) = 1.0. The 
green diamond denotes the averaged admittance at u) = 1.0. In contrast to the strongly disordered 
case Fig. IHl the circles concentrate in small region and the solid lines lie almost on the dashed line. 

random variables given by 



We call it the dielectric constant. It is obvious that the Cole-Cole plot of each sample is the 
Debye semi-circle since the form of Eq. fl23p is the same as that of the dielectric constant 
of the Debye relaxation. In addition, since |eAr(a;)| is less than unity, e]y{u) is finite for any 
set of the random jump rates Wi. 

In the analyses presented below, we use the following theorems of probability theory^ 
summarized in appendix [2 Let {xi} be a set of independent and identically distributed 
random variables. 

Theorem VI. 1 (Law of large numbers) When the expectation value (xi) exists, law of 
large numbers holds, i.e. 




(23) 



Prob 




as N 



oo. 



(24) 
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Theorem VI. 2 (Corollary of the extended central limit theorem) When the expec- 
tation value (xi) diverges, in order that the probability distribution ofY.f=iXi/aN — con- 
verges as N ^ oo it is necessary that the probability density p{xi) is of the form 

p{x) ~ x-^-^Lix) (25) 

for < a < 1. Here, L{x) is a slowly varying positive function which is precisely defined by 
Eq. \70^ . The normalization constant is chosen so that 

as N ^ oo. 

• When < a < 1 in Eq. probability distribution ofJ^^i/ciN converges to the 
stable distribution with the characteristic exponent a. 

• When a = 1 in Eq. probability distribution ofJ^^i/ciN — ^n converges to the 
stable distribution with the characteristic exponent a = 1. The centering constant 
diverges as N ^ oo^. 

By using these theorems, we discuss sample-dependence of the dielectric constant e{uj) = 
lim^r-^oo ^Afl"^)- The discussion is divided into two parts: the case where the expectation 
value (l/wi) is finite and the case where the expectation value diverges. 

When the expectation value (l/wi) is finite, it is shown from law of large numbers that 

1 1 

e(uj) = lim TT-, — ; — -- — = ; — ; — r- (27) 

Thus, there is no sample-dependence, i.e. the dielectric constant is self- averaging in the limit 
— i> oo. It is obvious that the Cole-Cole plot is the Debye semi-circle, which corresponds 

to the Cole-Cole plot of the dielectric constant for each sample. 

When the expectation value (l/wj) diverges (strong disorder), we discuss below the two 

cases of the value of a in Eq. (f^3j) : < a < 1 and a = 1. In order to see frequency 

dependence of e^u), the macroscopic limit (A^ oo) should be taken so that ujJ2'^/u!i/N 

is finite. Otherwise, e{u) becomes a constant or 1. 

• When < a < 1, from the extended central limit theorem, T = limjv_+oo J^'^/wi/aN is 
a random variable with a limiting probability distribution and hence it is finite with 
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FIG. 12: The frequency-dependence of the imaginary part of the disorder-averaged dielectric con- 
stant when the probabiUty distribution function of the jump rates is power-law p{wi) = aw"~^ (a = 
1/2). The two curves for sufficiently long chains = 100, 1000 lie on the same curve when plotted 
as a function of the scaled frequency ujN^^"'^^ {a = 1/2). 

probability one. Hence, the limit N oo should be taken so that uja^/N is finite; 

1 



e(uj) = lim 



= lim ■ ; ; — 

~ l + iuf ^^^^ 
It implies that the frequency-dependence is scaled with uoaN/N as shown in Fig. El 

Since T is a random variable, e(u;) is also a random variable and hence e(u;) is non-self- 
averaging. In addition, since it is known that the probability density p{x) of a stable 
distribution with characteristic exponent a behaves asymptotically as p{x) ~ x^"^^, 
the probability density p{T) behaves as p{T) ~ Hence, the expectation value 

(T) diverges. Since it implies that the coefficient of the first order of low uj expansion 
diverges, the disorder-averaged dielectric constant is a non-analytic function oi i uj 
which must behave as {t{uj)) — 1 ~ (i cD)^ (0 < ;U < 1). Thus, from the same argument 
employed in the discussion for the Cole-Cole plot of the admittance in Sec. |V| the 
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FIG. 13: The Cole-Cole plot of the disorder-averaged dielectric constant (solid line) when Wi are 
random jump rates with the power-law distribution Eq. 1)12(1 with a = 1/2 and the Debye semi- 
circle (dashed line). It is seen that the Cole-Cole plot of the disorder-averaged dielectric constant 
is inside the Cole-Cole plot for each sample, i.e. the Debye semi-circle. 

Cole-Cole plot of the disorder-averaged dielectric constant appears inside the Debye 
semi-circle which corresponds to the Cole-Cole plot of the dielectric constant for each 
sample (see Fig. IT!?jl. 

• When a = 1, from the extended central limit theorem, (wi) /dN — ^n is a random 
variable with a limiting probability distribution and hence it is finite with probability 
one. Furthermore, diverges as ^ oo as shown in appendix^ Hence, the limit 
N CO should be taken so that ujaNhN /N is finite; 

eiuo) = lim — ; — -- — 

1 



lim 

o—Loaj^h^ /N 
1 



1 + iuaNbN/N [(E(1M)/«7V - bN) /hn + 1] 



(29) 

Hence, e(u;) is self- averaging. 
Consequently, we conclude that only when < a < 1 the dielectric constant is non-self- 
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averaging. 

In order to test the foregoing discussion, we present an example of the above analysis 
where the random jump rates distributed by the power-law distribution Eq. (jl2|) when 
a = 1/2. From the extended central limit theorem, it is known that the distribution function 
of T is given by the stable distribution with a = 1/2. From Eq. (|8Up . the probability 
distribution g(T) of T behaves as g{T) ~ T^^/^. Hence, when u) ~ the disorder-averaged 
dielectric constant behave as 

(e(c<;)) - 1 ~ (30) 
On the other hand, the dielectric constant for each sample behaves as 

e{uj) — 1 ~ id;. (31) 

Thus, the low-frequency behavior is non-self-averaging. Absence of self-averaging is also 
seen in the Cole-Cole plot (Fig. IT^ . 

VII. HIGHER DIMENSIONAL DISORDERED MEDIA 

It is important to note that the system treated above is purely one-dimensional and it is 
still an open problem whether the admittance for higher dimensional media is self-averaging. 

In this section, we investigate properties of the admittance for higher dimensional disor- 
dered media in the same manner for one- dimensional media. At first, the statistics of the 
MFPT for higher dimensional media is analyzed analytically. With the results, it is shown 
that the complex admittance is non-self-averaging. After that, the theoretical results are 
confirmed numerically. 

For its feasibility, we investigate a site-disordered model called the random trap model, 
which is defined as a random walk model where the jump rate Wi j from site j to site i 
depends only on j. 

We assume that the higher dimensional medium is a x M square (or cubic) lattice. 
Here, A^ is the linear dimension of the medium in direction of the current and a number of 
grid points on the cross section perpendicular to the current direction is denoted as M. The 
periodic boundary condition on the side ends is assumed, i.e., the medium is a hyper-cylinder, 
and the oscillatory input current is introduced uniformly on the left end of the medium. An 
example of two-dimension is illustrated in Fig. El The probability distribution for the 
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FIG. 14: The schematic picture of the two-dimensional medium with an oscillatory input current. 

jump rates is assumed power-law distribution Eq. ()12|1 with the exponent a. Since we are 
interested in amorphous semiconductors, our analyses are limited to the case < a < 1. 

A. The analytic expression for the MFPT of higher dimensional media 

At first, we introduce a general analytical expression for the MFPT in terms of the 
backward equation. After that, by using the expression, the results for the one-dimensional 
MFPT and the extended central limit theorem, the probability distribution of the higher- 
dimensional MFPT is investigated. Since the analyses are general, results are applicable 
both to two-dimensional and three-dimensional media. 

The backward equation, i.e., the time-development equation for a averaged quantity, is 
derived in the following way. 

Let -Pji(t) be the probability that the random walker starting from the site i at time 
is found at the site j at time t. We assume that the master equation for the probability is 
written as 

Pi(t) = EV-lt). (32) 

Here, Pi is the vector whose j-th component is Pji with fixing the starting site i. Furthermore, 
if is a linear operator given in terms of the jump rates. 
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We consider the time-development equation of the average of the site-dependent quantity 
/j. The average gi is defined by 

^iW^E/j^JiW (33) 
j 

and hence the time-derivative is written as 

m = E/j^jiW = E/j^jk^kiW- (34) 

Taken the hmit t ^ 0, 

m = E/j^jk^ki = E/j^ji = E C^)y /j- (35) 

jk j j 

Let g(t) be the vector whose component is gi{t), then 

g(0) =* m. (36) 
By using the above result, g{t) is Maclaurin-expanded as 

sit) = E rfg^^no)^" = E ^,CHrit\ (37) 
Taken the time-derivative, 



n=0 n=0 ^■ 



OO in 



g(^) = E ^(*^)"fi"-' E ^ Hg{t). (38) 



n=i ^! ;;ro m 



The last equality gives the backward equation: 

g{t)^'Hgit). (39) 

As an example of the average gi{t), we consider the probability that the random walker 
staring at the site i on the end illuminated with the laser light is found in the medium at 
the time t{> 0). The survival probability ni{t) is defined as 

^i(^) = E^Ji(^)- (40) 
j 

Since the probability can be considered as a averaged quantity, ni{t) obeys the backward 
equation: 

h{t) Hn{t), (41) 
where n{t) is the vector whose component is ni{t). 
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In terms of the survival probability, the first passage time distribution (FPTD) F{t) is 
written as 

Fit) = -j^Y.^-.{t)/M, (42) 

i 

where M = X^i 1 and the sum is taken over the grid points on the end of the medium. Hence, 
the MFPT i is also written as 

roo 1 fOO rj 1 /-oo 1 

Here, ni(0) is the limit s — ^ of the Laplace transform n\{s). The Laplace-transform of the 
backward equation is written as 

sfi(s)-n(t = 0) =*i7fi(s). (44) 

Taken the limit s — > 0, 

*iffi(0) = -n(t = 0), (45) 

where n(t = 0) =* (1, 1, ■ ■ ■)• By solving Eq. (03), the MFPT is obtained as Eq. 

In the rest of this section, we concentrate on the random trap model. For the d- 
dimensional random trap model, the analytical expression of the MFPT can be easily ob- 
tained by Eq. ()45|1 in the following way. 

Each components in the same column of the matrix H are written only in terms of 
the jump rates of the same site. It means that each components in the same line of the 
transposed matrix *if are written only with the jump rates of the same site. Hence, divided 
the each lines of Eq. (jlSj) by the jump rates, we obtain 



*i^on(0) 



(46) 



V J 

Here, ^Hq is the transpose of the time- development operator when the jump rates of the all 
sites are equal to unity. From Eq. (|46|) . n(0) is a homogeneous linear function of the inverse 
of the jump rates. 

We denote the grid points in (i-dimensional medium with two indexes, n and m. One 
index n is the label of the cross sections normal to the current direction and the other index 
m is the d—1 dimensional vector to denote a grid point on the cross section. Let the cross 
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section at the left end, which is illuminated with laser light, be numbered and the cross 
section at the right end be numbered — 1. Each cross sections have M grid points. 

Since the boundary condition at the side ends is periodic, each grid points in the same 
cross section are identical geometrically and each jump rates Wnm in the same cross section 
n have the same contribution to the MFPT. Hence, the results obtained above means that 
the MFPT for the (i-dimensional medium, id is written as 

1 a 

m n=0 "^n-m 

= J_y(^ + ^ + ... + ^EzL\^ (47) 

where {a„} are constants independent of configuration of disorder. If all the jump rates in a 
cross section are equal, i.e., Wnm = Wnm' = Wn, the ci-dimensional MFPT is equal to the 
one-dimensional MFPT ii, which is given as 

N N-1 1 , , 

h = — + + ■ • • + . 48) 

Wo Wi Wn-1 

Comparison of Eq. ()47|) to Eq. PHj) shows that 

an = N -n. (49) 
Consequently, the MFPT for d-dimensional medium is given as 

M ^ \Wom Wim Wn -Ira) 

With respect to the MFPT, the d- dimensional medium is a bundle of one-dimensional chain 
and hence the MFPT is rewritten as 

~td = ^Y.t'r\ (51) 

where t^^^ is the MFPT for a one- dimensional medium defined as 

r.E^. (52) 

n=0 "^'im 
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B. The statistics of the MFPT for higher dimensional media 



In order to consider the probability distribution of td, the following postulate, which is 
expected to be true from the results of analyses for the one-dimensional case, is introduced: 
Let {un\ be a set of independent and identically distributed random variables. The tail of 
the probability density is assumed to be 

pu{Un) ~ tin~""\ (53) 

where < o; < 1. Then, we postulate that the random variable v defined as 

1 N-l 

jY.(N- n)u^ (54) 



n=0 

has a normalized probability distribution function and the tail obeys the power law. 

By making use of the postulate, we analyze the statistics of in two different ways and 
determine the asymptotic form of the probability distribution of t^- 

At first, we take the way where the sum in terms of n is done first. Prom our postulate, 
the random variable tra defined as 

= ^ (55) 

has the normalized probability distribution function and it obeys the power law: 

Pf(^m) - (56) 

where a' is a positive constant. Since id = M~^N^/"'^^Y,mim-i from the extended central 
limit theorem, the random variable x defined as 

^ = Ml/oc'-l]S[l/a+Jd (57) 

has a normalized probability distribution function and its tail asymptotically behaves as 

p,{x) ~ x-'^'-K (58) 

As the other way, we evaluate the sum in terms of m first. The extended central limit 
theorem shows that the random variable y„ defined as 

^""m^?^ (59) 
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has a normalized probability distribution function whose asymptotic form is 

PyiVn) ~ Vn-'-''. (60) 

Furthermore, since 1^ = M^/"~^ J2n=oi^ ~ ^)yni our postulate means that the random 
variable z defined as 

^ id (61) 



has a probability distribution function whose tail behaves as 

(62) 

Here, we compare the two ways of normalization Eqs. ()57p. ()61|). Since the two ways of 
normalization to obtain the normalized distribution function should be identical, 

a' = a. (63) 

Summarized our results, the MFPT td is normalized as 

td = M^/"-^ArV"+i;z (64) 

and the probability distribution function has the power-law tail as 

p^iz) z-'''\ (65) 

Furthermore, from Eq. ()63|). the tail of the probability distribution function of v defined by 
Eq. behaves as 

p,iv) ~ v-'^-K (66) 

Since the results are based on the postulate, we confirm numerically them. Figure IT5l 
shows that the probability distribution function of z converges to a unique distribution when 
N and M are increased. The power-law tail of the distribution is shown in Fig. ^1 It is 
clearly seen that the exponent is equal to —1.5 when a = 0.5. 



C. Non-self-averaging complex admittances for higher dimensional media 

By using the results obtained in the previous subsection, Eqs. (jMj) . we show in the 
same way as for the one-dimensional case that the complex admittance for higher dimensional 
media is also non-self-averaging when < a < 1. 
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FIG. 15: The probability distribution function of the normahzed MFPT t^A^^^/"^^M^^/"^^ for 
the random trap model where the probabihty distribution function of the random jump rates is 
power-law Eq. (|12|1 with exponent a = 1/2. The distribution is computed from 100000 samples 
of the length and width N,M = 2,5,10,100,300,1000. It is clearly seen that the probability 
distribution functions lie on a same curve when the length of lattice N is sufficiently large. 
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FIG. 16: The asymptotic behavior of the probability distribution function of the normalized MFPT 
computed from 1000000 samples. The solid line denotes the distribution function when N = M = 
100 and a = 1/2 and the dashed line represents the function proportional to z~^'^. One sees clearly 
that the distribution function behaves asymptotically as Pz{z) ~ a — 1 ~ 1.5 < 2). 
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The scaled frequency 



The scaled frequency 



FIG. 17: The frequency dependence of the real part (a) and the imaginary part (b) of the complex 
admittances for two-dimensional media when the exponent of the probabihty distribution of jump 
rates is a = 1/2 and the size of the media is = M = 10. The thin solid hues represent the 
admittances for each of 100 samples. The bold solid line represents the admittance averaged over 
20000 samples. Though the admittance of each sample is proportional to the scaled frequency uj, 
the averaged admittance is proportional to (/^ — 1/2)- 

The normalization of the MFPT, Eq. ()64p . means that the complex admittance for an 
infinitely large higher dimensional medium is given as 



Here, the infinite volume limit A^, M — > oo is taken with fixing uj at a finite value. Since 
the normalized MFPT 2; is a random number, the complex admittance is dependent on each 
samples even for the infinitely large medium. That is, the complex admittance is non-self- 
averaging. 

Furthermore, from Eq. it is shown that (z) = 00. It means non-analyticity of the 

complex admittance about u = and hence 



N,M^oo 




(67) 



in terms of the scaled frequency u defined as 



(68) 




(69) 



where /i is a constant such that < /i < 1 (see Fig. [T7j). 
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FIG. 18: (a) The Cole-Cole plot for 4x4 disordered square lattice. The probability distribution 
for the jump rates is power-law distribution Eq. (|12() when a = 1/2. The solid lines denote 
admittances for 100 samples and the dashed line denotes the admittance disorder-averaged over 
5000 samples, (b) The Cole-Cole plot for 10 x 10 disordered square lattice, (c) The Cole-Cole plot 
for 30 X 30 disordered square lattice. These numerical results show that the sample-dependence of 
the admittance seems to exist even when the size of the lattice is large and hence it means that 
the admittance is non-self-averaging. 

The results obtained above are all concerned with the low frequency behavior. We in- 
vestigate numerically the non-self-averaging characteristics in the all frequency range with 
the Cole-Cole plot. Figure UHl shows the Cole-Cole plots of the admittances for one hundred 
realizations of 4 x 4, 10 x 10, 30 x 30 square lattices. These numerical results show that the 
sample-dependence of the admittance seems to exist even when both the length of the lattice 
and the width are large and hence it implies that the admittance is non-self-averaging even 
in higher dimensional cases. 
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VIII. DISCUSSION AND CONCLUSIONS 



In conclusion, we presented a generic stochastic dynamical model in one-dimensional 
medium of the BPM and showed that the admittance has important information about 
transport in disordered media, e.g. the FPTD. Thus, the BPM will be a powerful technique 
to investigate the dynamical process in random media. We also showed that the admittance 
is non-self-averaging when the jump rates are random variables with the power-law distribu- 
tion as p{wi) ~ (0 < a < 1). It is important to note that this power-law distribution 
is realized in amorphous semiconductors. 

On the other hand, when 1 < a, the low frequency behavior of 0{uj) has no sample 
dependence as shown in Eqs. (II 7|) . (fT5|) . However, with respect to the terms of 0(^7^), 
the second moment of the inverse of the jump rate {l/w"^), which diverges when 1 < a < 
2, is concerned. It implies that some anomalous statistical properties may be observed. 
Properties for the full frequency range were analyzed with the Cole-Cole plots. In the 
marginal case when a = 1, the non-self-averaging property appears as shown in Fig. ITnl 
However, Fig. 1111 shows that such effects from terms of 0(u;^) are small when 1 < a. 

When the complex admittance is non-self-averaging, the behavior of the averaged admit- 
tance over infinite number of samples does not coincide with that of the admittance observed 
in experiments. In order to see the non-self-averaging behavior described above one needs 
to work in one- dimensional system or the case where the translational invariance is violated 
in only one direction. 

Furthermore, in order to investigate possibility that the non- self- averaging properties of 
the admittance are observed in real experiments on higher dimensional media, we performed 
theoretical and numerical analyses of the admittance on higher dimensional disordered me- 
dia. The analyses on anomalous statistics of the MFPT show that the non-self-averaging 
characteristics are observed in the low-frequency behavior of the complex admittance of 
higher- dimensional media. The numerical results confirm the conclusion that the admittance 
is non-self-averaging. Hence, it suggests possibility that the non- self- averaging properties of 
the admittance of the BPM may be observed in the real experiments on higher dimensional 
media. 

On the other hand, non-self-averaging behavior, i.e. strong sample-dependence may be 
observed in other experiments. For example, it is known that anomalous system-size depen- 
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dence of the mobility, fi ~ L^~^/°', is observed in the time-of-flight experiment of amorphous 
semiconductors^. The anomalous system-size dependence is due to the fact that the small- 
est value of random jump rates depends on the length of a chaini^. Since the smallest 
value dominates the mobility, the mobility is also a random variable. Thus, the anomalous 
system-size dependence of the mobility implies that the mobility is not self-averaging. 
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LAW OF LARGE NUMBERS AND THE EXTENDED CENTRAL LIMIT THE- 
OREM 

In this appendix, we summarize the notions of probability theory, partly rewritten for 
our convenience, used to analyze the complex admittance and the "dielectric constant" in 
Sees. El and ED After that, we show that Eq. (fT!?jl diverges as N oo. 

1. The extended central limit theorem 

Here, we summarize theorems on probability distribution of normalized sum of random 
variablesti^. 

The following law of large numbers tells that a mean converges to an expectation value 
even when a variance does not exist. In this section, S'„ denotes the sum J27=i -^i- 

Theorem .1 (Strong law of large numbers) Let the Xk be independent random vari- 
ables with a common distribution F. If they have an expectation fi then Sn/n —>■ /i with 
probability one. 

In order to describe the extended central limit theorem, some concepts are defined. 
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Definition .1 (Domain of attraction) The distribution F of the independent random 
variables Xk belongs to the domain of attraction of a distribution R if there exist the norming 
constant a„ > and the centering constant bn such that the distribution of a~^Sn — &n tends 
to R. 

Definition .2 (Slow variation) A positive function L defined on (0, oo) varies slowly at 
infinity if 

for each t > as s oo. 

It implies that the leading term of L{st) is L{s) when s is large. Constant functions and 
log(a;) are examples. 

In addition, the function F{x) is assumed to be related to the probability density p{x) as 



F{x) ^ / dyp{y). (71) 

J — oo 

Using definitions presented above, the extended central limit theorem is stated as follows: 

Theorem .2 (The extended central limit theorem) In order that F belongs to some 
domain of attraction it is necessary that F is of the form 

1- F{x) +F{-x) x-''L{x) (72) 

for some 0<q:<2 as x ^ oo. 

• If a = 2 then F belongs to the domain of attraction of the normal distribution. 

• If a < 2 and 

l-^M p. , J'-^l (73) 



1- F{x) + F{-x) 1- F{x) + F{-x) 

as X ^ oo then F belongs to the domain of attraction of the stable distribution with the 
characteristic exponent a. If Eq. fails then F belongs to no domain of attraction. 

Here, Eq. ()73p only means that the limits exist. If positive random variables are considered, 
p = 1 and g = since F{—x) = 0. The limits p and q give parameters of a stable distribution. 
In addition, the centering parameter 6„ is given as followings; 

• If a > 1 then 6„ is given by the expectation value. 
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If a = 1 then 6„ is given by 

' X 



oo 

bn = n I dxsm[ — ) p{x). (74) 

(In , 



oo 



If a < 1 then 6„ = 0. 



The norming constant a„ is chosen so that 



77/ 

L{an) ^ C, (75) 



a 

"n 



where C > is a constant. 

This theorem implies that "when the disorder is strong and the probability distribution 
function of the normalized sum Sn/^n — bn converges in the limit n oo the limiting 
probability distribution is the stable distribution with the characteristic exponent < a < 1 
given by Eq. f7^ . The inequality < a < 1 is due to divergence of the expectation value 
(the strong disorder). 

For our convenience, we rewrite the necessary condition Eq. ()72|) for the convergence of 
the probability distribution of S'„/a„ — 6„ in terms of the probability density p(x). Since 
p(x) = F\x) = dF{x)/dx, the asymptotic behavior of the probabihty density is given by 

pix) + pi-x) ~ a:-"-iL(x)[l - -^^]. (76) 

The second term on the right hand side is evaluated from Eq. (f70|l . By differentiation with 
respect to t of Eq. (fTUI) and setting t = 1, it is shown that 

hm ^ = 0. (77) 

It implies that 

p{x) + p{-x) ~ a;-"-^L(x) (78) 

as X oo. 

For the asymptotic behavior of the stable distribution, the following theorem is known. 

Theorem .3 (Asymptotic behavior of a stable distribution) If G is a stable distri- 
bution obtained as the limiting distribution of the extended central limit theorem, then as 
X — > oo 

2 — 0' 2 — (T 

x"[l - G{x)] ^ Cp -, x'^Gi-x) Gq -, (79) 

a a 

where the parameters G,p,q are given by Eqs. and |73D . 
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It implies that the asymptotic behavior of the probabihty density g{x) = G'{x) is given by 

g{x) - Cp{2 - g{-x) ~ Cq{2 - a)x-°-\ (80) 

Hence, a stable distribution with the characteristic exponent a has absolute moments of all 
orders < a and absolute moments of all orders > a do not exist. 

Next, we present the relation holding when < a < 1. It tells about the contribution 
of the maximal term to the sum. Let Xi,X2, ... be independent random variables with 
the common distribution F satisfying Eq. ()72|) with < a < 1, i.e., belonging to the 
domain of attraction of the stable distribution with the characteristic exponent a. Put 
Mn = max(Xi, . . . , X„). Then, 

as n ^ oo. It implies that the maximal term is of the same order as the sum with probability 
one when < a < 1. 

2. The centering constant 6„ when a = 1 

Here, we show that the centering constant 6„ diverges for the positive random variables 
X satisfying Eq. (f?^ with a = 1; 

1- F{x) x-^L{x). (82) 

Let Xq > denote the lower limit of the distribution and the probability density is given 
by p{x) = F'{x). From Eq. (j75p . the norming constant a„ is chosen so that 

T) 

-L{a^) ^ C, (83) 

where C > is a constant. The asymptotic behavior of the probability density p{x) is 
derived from Eq. (fTHj) : 

p{x) ~ x~^L{x). (84) 
We evaluate 6„ by dividing the region of integration of Eq. (f7i|) ; 

bn = n dxsin( — j p{x) + n dxsm( — j p{x) (85) 

where x„ is a constant satisfying 

Xn "C fln, lim Xn = OO. (86) 
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Since x„ <C a„, the first term on the right hand side of Eq. (jH^ is given by 

the first term ~ — / dxxp{x) > 0. (87) 
From Eq. (jHl}, the leading term of the second term is evaluated as 



oo 



the second term ~ n / dxsin ( — ] x L(x) 



X 



/ sin (x) x~^L(a„x), (88) 

la n. 



n 
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where a numerical constant is ignored. Since L is a function slowly varying, the leading 
term is given by 



oo 

-2 



the second term ~ — L{an) I dxsm{x)x' 

dn •Ixn/an 
Tl 

L(a„)log(a„). (89) 

On 

From Eq. (jHHj) . it is shown that 

the second term ~ log(a„). (90) 

Consequently, since the first term is positive and the second term becomes positive infinity 
in the limit n — > oo, the centering constant diverges in the limit n — > oo. 
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